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Abstract 

We formulate in terms of the quantum inverse scattering method the algebraic Bethe 
ansatz solution of the one-dimensional Hubbard model. The method developed is based 
on a new set of commutation relations which encodes a hidden symmetry of 6-vertex 
type. 



In 1968 the exact solution of the one- dimensional Hubbard model by the coordinate Bethe 
ansatz was presented by Lieb and Wu |2[ . It took some years until Shastry found the 
many conserved charges|| and also the two-dimensional classical vertex model^, |5j whose 
transfer-matrix commutes with the Hubbard Hamiltonian. The .R-matrix responsible for the 
integrability( "infinite number of conserved charges") was then explicitly exhibited [|], [|. Some 
time later, Wadati et al were able to verify such results by using a quite different and 
interesting approach. Afterwards, Bariev developed a variant of the coordinate Bethe 
ansatz to study Shastry's vertex model, though on the basis of a diagonal-to-diagonal transfer 
matrix approach |J. More recently, some progress has been made concerning the Yangian 



symmetry of the Hubbard model M and also on its "free-fermion" Yang-Baxter structure [10 



Some further discussion about Hubbard's invariant can also be found in the literature 11 



However, certain important properties underlying such "integrable" program still needs to 
be understood. This is justified, for example, by the early attempt of Shastry || in conjecturing 
the eigenvalues of the row-to-row transfer matrix of the "covering" vertex model. An important 
step towards closing this program is certainly the formulation of the Bethe states of the one- 
dimensional Hubbard model by means of the quantum inverse scattering approach |13| . Unlike 
the standard Bethe ansatz, this method is based on first principle algebraic rules and definitely 
brings new insight on the mathematical structure of integrable systems. The solution of the 
Hubbard model by the quantum inverse scattering method is, in fact, a long-standing problem 
in the field of exactly solved models. In this letter we show how this more unified approach 
of Bethe ansatz technique can be established for the one-dimensional Hubbard model. In the 
course of our formulation we had to overcome a major difficulty: the non-additive property 
of the Hubbard R matrix. We have found the fundamental commutation rules between the 
creation and annihilation operators present in the embedding vertex model. It turns out 
that the eigenvectors, eigenvalues and the Bethe ansatz equations follow as a consequence of 
systematic algebraic manipulation of such commutation rules. A hidden symmetry of 6- vertex 
type, important for integrability, is noted. We think that our results should also be of relevance 
for future developments of the physical properties of the one-dimensional Hubbard model. One 
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possibility should be the application of our formulation in the context of Korepin et al. method 
|T5|j of computing correlation functions. 

The Hamiltonian of the Hubbard model on a one-dimensional lattice of length L is written 
in the form 

L L 

^ = "E E (ctaQ+v + h.c) + niA n i,l C 1 ) 

i=l cr=T,i i=l 

where c; )0 -(nj jCr ) are Fermi (number) operators with spin a on site i, and U is the Hubbard 
coupling. Fundamental to the integrability of the Hubbard model is the fact that Hamilto- 
nian (1) commutes with certain one-parameter family of transfer matrix T. In analogy with 
integrable systems in classical mechanics, T is the generator of the many conserved charges, 
and Hamiltonian (1) is one of those currents. Thus, the analysis of the physical properties of 
the transfer matrix will certainly provide a deeper understanding of the Hubbard model. The 
appropriate two-dimensional classical statistical system exhibiting such properties was found 
by Shastry ||, [5]. The model is constituted of two coupled 6- vertex model satisfying the free- 
fermion condition. The vertex model is parametrized in terms of three functions a(X),b(X) and 
h(X) which are constrained by Hubbard's coupling as 

sinh[2h{X)\ = -a{X)b{X) (2) 

where functions a(A) and 6(A) are the non-trivial free-fermion Boltzmann weights. This gives 
us an i?-matrix R(X, fi) consisting of ten distinct Boltzmann weights. Here we denote them by 
ai(X,fi),i = 1, • • • , 10. In Appendix A we present the structure of the i?-matrix, the explicit 
expressions and some usefull identities for the weights «j(A, n) @, |J. In general, the transfer 
matrix T is obtained as a trace of an auxiliary monodromy operator, T = TtqT . The space G 
is a " ghost " variable, corresponding to a horizontal arrow in the classical vertex model. Its 
dimension corresponds to the four possible states of the Hubbard model on a given site. As 
we shall see below, it is convenient to write the associated monodromy matrix T(A) as 

/B(X) B(X) F(X) \ 
T(A)= (7(A) A(X) B*(X) (3) 

\C(X) C*(X) D(X)J 



where B(X) (B*(X)) and C(X) (C*(A)) are two component vectors with dimensions 1 x 2(2 x 1) 
and 2 x 1(1 x 2), respectively. The operator A(X) is a 2 x 2 matrix and the other remaining 
operators are scalars. The integrability condition is based on the Yang-Baxter algebra, namely 

R(X, n)T(X) <£> T(ji) = T(jm) ® T{X)R(X, //) (4) 



where the symbol ® stands for the Grassmann direct product ]T2|. Such definition takes into 



account the extra signs appearing when fermionic states ( spin up and down ) are permuted 
||. One consequence of Shastry's Boltzmann weights is that the monodromy matrix has a 
triangular form when acting on the standard ferromagnetic pseudovacuum |0). More precisely, 
we find the following diagonal properties 

B(X) |0> = l^e 2h ^] L |0> , D(X) |0> = [^|e 2 ^] L |0> , A aa {X) |0) = |0) , a = 1, 2 (5) 

as well as the annihilation identities 

C(X) |0) = (7(A) |0> = C*(X) |0) = 0, A ab (X) |0) = 0(a ^ b = 1, 2) (6) 

This suggests that operators B(X), B*(X) and F(X) act as creator fields on the ferromagnetic 
reference state |0). We notice, however, that the operators B(X) and B*(X) do not mix under 
the integrability condition (4). Therefore, in the construction of the eigenvectors it will be 
enough to look only for combinations between the fields B(X) and F(X). The one-particle 
state |$i(Ai)) is made by the linear combination 

\^ 1 (X 1 )) = ^(Ai).^ |0> = B a (X x )T a |0> (7) 

where T a is the component of a constant vector JF with dimension (2x1). The two-particle state 
|$2(Ai, A2)) depends both of operators B(X) and F(X). This happens because the commutation 
rule between two fields of type B(X) generates the scalar operator F(X). This is a constrain 
imposed by the integrability condition (4), which reads 

B(X)®B(V) = a 1)2 (A,/x)[B(/x) ® B(X)}.r(X,fi) - ia 10t7 (X,fi){F{X)B(ji) - F(fi)B(X)}£ (8) 
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where we define a a ,&(A, aO = ot a (X, /i)/«b(A, //). The vector £ and the matrix r(A,/i) have the 
following structures 

/l 0\ 

a(A,/x) 6(A,/i) 



e= (0 1 -1 0), f(A,/i) 



6(A,/i) a(A,/i) 
\0 1/ 



(9) 



where 



(10) 



6(A, n) = a 8j i(A, //)a 9)7 (A, //), a(A, /i) + 6(A, //) = 1 

Remarkable enough we have found that f-matrix (9) is in fact factorizable. Moreover, when 
properly parametrized, it has the same structure of that appearing in the isotropic 6-vertex 
model. We stress that such hidden symmetry is crucial in our algebraic construction and plays 
a fundamental role on the exact solution of the Hubbard model. In our opinion, this is the 
" nice " algebraic explanation for the fact that the bare two-body scattering of the Hubbard 
Hamiltonian appears in the 6-vertex form [l], |^| . This result can be established by performing 
the following change of variables 



A 



o(A) 



,2ft(A) 



6(A) 



2h(X) 



u 

~2 



b(X) a(A) 

By using the Boltzmann weights [|], || ( see Appendix A ) in equation (10) and by con- 
sidering the new variables defined in (11), we are able to rewrite functions a(A,/i), b(X,fi) 
as 



a(A,/i) 



U 



fi-x + u' 



b(X,jl) 



pL — X 

jl-X + U 



(12) 



which are precisely the non-trivial Boltzmann weights of the isotropic 6-vertex model [13|, 
|i~4l , Taking into account our considerations above, it is not difficult to check that the 
two-particle state is given by 



|$ 2 (Ai, A 2 )> = {B(Ai) ® B(X 2 ) + iai 0l7 (Ai, X 2 )F{X 1 )(^ ® %)B(X 2 )}.f |0) 



(13) 



where $o is the unitary constant. In fact, we have checked that all unwanted terms generated by 
the eigenvalue problem can be canceled out through a unique Bethe ansatz equation. Moreover, 



at least at this level, the physical meaning of our construction is the following. While each 
component of the field B(X) creates an electron with spin up or down, the operator F(X) 
is responsible for the double occupancy on a given site of the lattice. In general, the n- 
particle state can be constructed by induction and we have verified that it satisfies the following 
recurrence relation 

|$n(Al,-",An)) =Q n (\ 1 ,---,\ n ).?\0) (14) 



where 



$n(Ai, • • • , A n ) = B(X 1 ) <g> $ n -i(A 2 , • • • , A„) + «aio,7(Ai, Xj) J[ ia 2 ,9(X k , Xj) 

j=2 k=2,k^j 



X 



|*(8) F(Ai)$ n _ 2 (A 2 , • • • , Aj_i, Aj+i, • • • , X n )B(Xj)\ J| ai, 2 (A fc , Xj)f k ,k+i(Xk, Xj) (15) 

k=2 

Let us now turn to the diagonalization problem. The associated transfer matrix is obtained 
as a graded trace of the monodromy matrix T(A). The graded structure takes into account 
the fermionic degrees of freedom, and on the diagonal of T{X) only A aa (X) contributes with a 
non-null Grassmann parity. Hence, the eigenvalue problem becomes 

[B(X) - £ A aa (X) + D(X)] |$ n (Ai, • • • , A n )) = A(A, {A,}) |$ n (Ai, • • • , A„)) (16) 

a=l 

In order to solve (16) we need the commutation rules between the diagonal and the creation 
operators. This is similar to solving a problem of quantum mechanics on the Fock space, 
analogously to the role of the Heisenberg algebra on the solution of the harmonic oscillator. In 
our case, the necessary commutation relations can be obtained by an appropriate manipulation 
of integrability condition (4). The procedure is rather cumbersome, and here we only list some 
of them which are fundamental for further discussion. They are given by 

A(X)®B(fi) = -ia h9 (X,n)[B(fi) ® A(\)].r(\, n) +ia 5fi (X, fi)B(X) ® 
{-iai 0)7 (A, n) [B*(X)B(n) + m 5 , 9 (A, fj)F(X)C(jj) - m 2 , 9 (A, /i)F(/i)C(A)] } ® £ (17) 

B(X)B(fj) = ia 2i9 {v, X)B(fj)B(X) - m 5 , 9 (/i, X)B(X)B(fx), (18) 
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D(X)B(ii) = -ia Si7 {\,v)B(ii)D(X) + a 5 j{\,ii)F{u)C*{X) 

-a 4 , 7 (A,/x)F(A)C*(/x) - ia 10>7 (\,n)(.[B*(\) ® i(/i)] (19) 

The eigenvalue A(A, {Aj}) can be calculated by keeping the terms proportional to the eigen- 
state |$ n (Ai, ■ ■ • , A„)). For example, by using several times the first terms of the commutation 
relations (17, 18, 19) we find the following structure for eigenvalue A(A, {Aj}) 

A(A, {A,}) = [^e 2 ^ ft fcMAi, A) + M e ^] L f[ -ia s , 7 (X, A,) 
& (A) t=i a ( A ) i=i 



- n -*ai, 9 (A, Ai)A (1) (A, {A,}, {^}) (20) 

8=1 

where A^^(A, {A.j}) is the eigenvalue of an auxiliary inhomogeneous problem related to the 
hidden 6-vertex symmetry we have mentioned before. More precisely, such auxiliary problem 
is defined by 

«(A, A 1 )f b d - 2 (A, A 2 ) • • -ffcrU \n)^- ai = A«(A, {A,}, { N })F b ^ (21) 

Fortunately such additional eigenvalue problem can be solved using the well known results 
of Faddeev et al |I3fl . New parameters {/ij} are then introduced in order to perform the 



diagonalization problem (21) . Here we just have to adapt their algebraic results in order to 
consider the 6-vertex problem on an irregular lattice. Considering that this later eigenvalue 



problem has appeared in many different contexts in the literature [OL n, 15], we just present 
our final results. First it is convenient to generalize a bit Shastry's parametrization []5J by 
introducing the new functions z±(x) as 

z _{x) = jMc^, z + (x) = Mg2fc(») (22 ) 
b(x) a(x) 

In terms of functions z±(x) and the variables {/!•/} introduced in (11), we find that the 
eigenvalue (20) (modulo overall constant)can be written as 

" 6(A) 1 + z-(Aj)/z + (A) " s_(A) - l/z_(A) - + U/2 
l\ o(A) 1 - z.{\)/x.{\) /i z_(A) - 1/*_(A) - ft - U/2 
_ » 6(A) l + ^-(AQ^(A) ™ l/z + (A)-^(A)-^--[//2 
^ g(A) 1 - MA*K(A) /i l/*+(A) - z+(A) - ft + t//2 l J 



Analogously, in order to cancel the unwanted terms, it is possible to show that the nested 
Bethe ansatz equations constraining the numbers {Ai},{/L,} are then given by 

" s_(A fc ) - l/z4X k ) - g, - U/2 _ " h- gj + g 

fe i = i 1 ^(A fe )-l/^(A fe )-^ + C//2 fAiLi-h-V 1 J 

To check the consistency of our results (23-25) one has to verify that A(A, {z±(Xi)}, {/%}) is 
free of poles for finite values of A. In fact, the null residue condition on both direct ( Z- (A) ) and 
crossed (z+(X) ) poles lead us to the Bethe conditions (24,25). A possible physical application 
of the eigenvalue result (23) is probably concerned with the finite temperature properties of 
the one-dimensional Hubbard model |jl~6| . This is connected to the recent developments of 
new powerful methods to deal with finite size effects in integrable models (T7|, P|. These 
techniques depend much on the diagonalization of the quantum transfer matrix (rather the 
one-dimensional Hamiltonian) , a problem which we managed to solve in this letter. Lastly, it 
is also possible to rewrite the nested Bethe ansatz equations (24,25) in terms of the original 
form presented by Lieb and Wu ||]. In this case, one just needs to change jlj — > 2ip,j and relate 
the variables A& with the lattice momenta p k by the relation z_(A^) = e %Vh . 

We would like to conclude this letter with the following comments. The eigenvalue (23) 
is almost the one conjectured by Shastry in ref. ||. They differ by important phase factors, 
which are not easily obtained by using only phenomenological arguments. Our result (23-25) is 
connected with periodic boundary conditions, while that conjectured by Shastry is related to 
a rather peculiar (sector dependent) toroidal boundary conditions. The method we presented 
here is easily extended for a more general inhomogeneous model [[|, ICfl . We expect that the 
only change in the Bethe ansatz equations (24,25) will be on the terms proportional to the 



power of L. We plan to discuss these results in a more detailed version of this letter ||19|| . 
Finally, some extra remarks are now in order. It is possible to show, from the commutation 
rules between the "dual" field B* (A) and F(X), that a second SU(2) 6- vertex hidden symmetry 
is also present Jl9|]. Thus, the two 6- vertex structure are tied up by the same field F(X). This 
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resembles much the constrain leading to the SO (4) symmetry of the Hubbard chain [ 2(| . This 
is known to be of relevance for the Bethe ansatz completeness ||21| , for the classification of 



the elementary excitations p2L and can also play an important role in the computation of 



correlation functions [15 
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Appendix A : 

The structure of the .R-matrix 0, || || is 
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1 a 2 
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"2/ 



where the weights aj(A,/i) ( normalized by a$(\, jj) ) are given by 

= {e [w ~ h(A)1 a(A)a(/i) + e- [hM ^ (A)] 6(A)%)}a 5 (A,/i) 
i(A,//) = {e- [hM - /l(A)] a(A)a( / u) + e [hM - h(A)] 6(A)%)}« 5 (A,/i) 

a 3 (A,/i) : 



«2l 



[MM)+MA)] Q ( A ) & ( /i ) + e -[ fe M+MA)] & ( A ) a ( /i ) J- C osh[/i(/i) - h(X)} \ 
a(A)6(A) + aQfibQi) \cosh[/i(/i) + h(X)] J 



(A.l) 



(A.2) 
(A.3) 
(A.4) 
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. _ e-^ +h ^a(X)b(fi) + e^ +h ^H(X)a(ij) ( cosh(/i(/x) - /i(A)) 1 
MX ^ ] ~ a(A)6(A) + a(/x)%) \ cosh(/i(/x) + /x(A)) J (A ' 5) 

(A.6) 

f -e-[ fe W+MA)] a(A)&(/i) + e [^)+MA)] 6(A)a(/x) | ^^ coshK/x)-/^)] 

a7(A '^ = 1 a(A)6(A) + afr)6fr) 1 [b ^ ] - h ^ cosh^ + MA)] " 5 ^^ 

(A.7) 

a 8 (A, /i) = {e [/lM - h(A)] a(A)6( / u) - e~ [/l{ ' l) - /l{A)1 6(A)a( / u)} a 5 (A, /i) (A.8) 

a 9 (A, /i) = {-e- [hM - hW] a(X)b(n) + e [h ^- h{x)] b(X)a(fi)} a 5 (A, /x) (A.9) 

/x x ^ 2 (/w)-& 2 (A) fcosh[/i(/x) -/i(A)]l , , 

aio(A ' ^ = a(X)b(X) + a(^) { cosh^xx) + MA)] J MA ' ^ ( ° } 

We remark that we have used the original Shastry's Boltzmann weights [Q, HJ together with 
the grading modifications of Wadati et al ||. Moreover, the 6- vertex parameters a(A) and 6(A) 
satistfy the free-fermion condition a(A) 2 + b(X) 2 = 1. We also list some important identities 
between the Boltzmann weights || 

a 3 (X, /x) = «i(A, /i) + a e (X, /x); a 4 (A, /i) + a 7 (A, /x) = a 2 (A, /x) (A. 11) 

a 2 (A, /i)«i(A, /i) — a$(X, /i)as(A, /x) = a 4 (A, /x)«3(A, /x) — «i (A, /x) = aj^A, /x) (A. 12) 
a 2 (A, /x)a 3 (A, /x) + a 4 (A, /x)ai(A, /x) = 2a;?(A, /x) (A. 13) 
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